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Abstract. Stieltjes' work on continued fractions and the orthogonal polynomials 
related to continued fraction expansions is summarized and an attempt is made to 
describe the influence of Stieltjes' ideas and work in research done after his death, 
with an emphasis on the theory of orthogonal polynomials. 
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Introduction 

The memoir Recherches sur les fractions continues, published posthumously in 
the Annales de la Faculte des Sciences de Toulouse — a journal of which Stieltjes 
was one of the first editors — and a great number of other papers by Stieltjes con- 
tain a wealth of material that still has a great impact on contemporary research, 
especially on the theory of orthogonal polynomials. The general theory of orthog- 
onal polynomials really started with the investigations of Chebyshev and Stieltjes. 
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The impact of the work of Chebyshev and his student Markov has already been 
described by Krein [53] . Here we give an attempt to discuss some of Stieltjes' con- 
tributions and the impact on later work. Orthogonal polynomials offer a variety of 
results and applications. The bibliography [91] up to 1940 consists of 1952 papers 
by 643 authors. Even now interest in orthogonal polynomials is enormous. One 
of the reasons is that orthogonal polynomials seem to appear in a great variety of 
applications. Their use in the numerical approximation of integrals was already 
pointed out by Gauss and further extended by Christoffel [19] [30] and Stieltjes 
[95] . The Pade table [73] for the approximation of a function by rational functions 
is very closely related to continued fractions and Stieltjes' work may be considered 
as one of the first proofs of convergence in the Pade table [74]. In 1954 Lederman 
and Reuter [56] and in 1957 Karlin and McGregor [49] showed that the transition 
probabilities in a birth and death process could be expressed by means of a Stieltjes 
integral of orthogonal polynomials. Even in pure mathematics there seems to be a 
natural framework where orthogonal polynomials come into play: representations 
for certain Lie groups very often are in terms of special functions, in particular or- 
thogonal polynomials (see e.g. Vilenkin [117]). Recently this has also been observed 
for quantum groups [52]. Discrete orthogonal polynomials have useful applications 
in the design of association schemes and the proof of nonexistence of perfect codes 
and orthogonal polynomials on the unit circle have a close connection with digi- 
tal signal processing. The proceedings of the NATO Advanced Study Institute on 
"Orthogonal Polynomials and their Applications" (Columbus, Ohio 1989) [69] gives 
excellent contributions to each of these aspects of orthogonal polynomials and is 
strongly recommended. 

Stieltjes' work has already been discussed by Cosserat [22] shortly after Stieltjes' 
death in 1894. In these notes we will try to estimate the value of the investigations 
by Stieltjes a century later. Let me also mention Brezinski's book on the history 
of continued fractions [10, Chapter 5, Section 5.2.4 on pp. 224-235] where Stieltjes' 
work on continued fractions is shown in its historic context. 



1. Stieltjes Continued Fraction 
The object of his main work [105] is the study of the continued fraction 



(1.1) 



CiZ + 



C2 + 



c 3 z H h 



C2n + 



1 



which is nowadays known as a Stieltjes continued fraction or S'-fraction. Stieltjes 
only considers the case where Ck > (k = 1, 2, . . . ). In general an S'-fraction is any 
continued fraction of the form (1.1) in which all are different from zero or any 
continued fraction which can be obtained from it by an equivalence transformation 
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contraction to a J-fraction 



(1.2) C - 
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with a 2 , = 1/ci, 60 = — l/(ciC2) and 
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a n — 2 7 bn — ■> k — 1,2, 



The positivity of the Ck, as imposed by Stieltjes, clearly puts some constraints on 
the coefficients (ik,bk e.g., a\ > and bk < 0. A J-fraction can be regarded as 
being generated by the sequence of transformations 

1 a 2 

t (w) = —, t k (w) = z-b k -i -, fc = l,2, ... 

u; «; 

The superposition to(^i(' ' ' (^n(w)) • • • )) for w = 00 is then the n-th approximant 
or n-th convergent of the fraction (1.2). This is a rational function of the variable 
z and we have ^ 

f (fl(-"(fn(oo))...)) ' / '"- |i ' J 



ai p n {z) 

where both the denominator polynomials p n (z) (n = 0,1,2,...) and numerator 
polynomials p < ^l 1 (z) (n = 0,1,2,...) are solutions of the three-term recurrence 
relation 

(1.3) zr n (z) 

— a n+l r n+l 

with initial condition 

p_i(z) = 0, p (z) = 1, p { l\{z) = 0, p$\z) = 1. 

The convergents of the S- fraction are such that the 2n-th convergent of (1.1) is 
equal to the n-th convergent of (1.2). If the denominator p n (z) vanishes for at most 

a finite number of integers n and if lim^^oo p ( ^}_ 1 (z)/p n (z) = f(z) exists, then the 
J-fraction converges to f(z). Stieltjes gave a general theory of S'-fractions (and 
consequently of J-fractions) with cj~ > (k = 1, 2, . . . ), dealing with questions of 
convergence and he showed a close connection with asymptotic series in terms of a 
given sequence of moments (see also the next section). 

One of the most important facts in the theory is that the denominators p n (—x) 
(n = 0, 1, 2, . . .) form a sequence of orthonormal polynomials on [0, 00) i.e., there 
is a positive measure fx on [0, 00) such that 



poo 

/ p n (-x)p m (-x) du(x) = S m . n . 
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The support of the measure fi is in [0, oo) precisely because Stieltjes assumes the 
coefficients Ck of the S- fraction (1.1) to be positive. Stieltjes showed that such 
orthogonal polynomials have zeros with interesting properties. He proved that all 
the zeros of p n (— x) are real, positive and simple; moreover the zeros of p n {— x) 
interlace with the zeros of p n -\{— x) but also with the zeros of P^-i{~ x). The 

latter property shows that the convergent p^-i( z ) I 'Pn( z ) is a rational function with 
n real and negative poles and positive residues. These properties are now quite 
classical and of great use for numerical quadrature. The property of orthogonality 
is crucial in these considerations (but Stieltjes never uses this terminology). A 
famous and very important result in the theory of orthogonal polynomials on the 
real line is the following result: 

Theorem. Suppose a system of polynomials satisfies a three-term recurrence re- 
lation of the form (1.3) with ak+i > and bk G K (k = 0,1,2,...) and initial 
conditions r-i(z) = and r (z) = 1, then these polynomials are orthonormal in 
L 2 ((i) for some positive measure fx on the real line. 

This theorem is usually called Favard's theorem [25] but it is basically already 
in Stieltjes' memoir [105, §11] for the case of J- fractions obtained from contracting 
an S'-fraction with positive coefficients: he shows that there is a positive linear 
functional S such that S(r m r n ) = whenever m ^ n. The only thing that Stieltjes 
was missing was the Riesz representation theorem which would enable one to express 
the linear functional S as a Stieltjes integral. 

Hilbert's work on quadratic forms in infinitely many variables was much inspired 
by Stieltjes' work on continued fractions [43, p. 109]: "Die Anwendungen der The- 
orie der quadratischen Formen mit unendlich vielen Variabeln sind nicht auf die 
Integralgleichungen beschrdnkt: es bietet sich nicht minder eine Beriihrung dieser 
Theorie mit der schdnen Theorie der Kettenbruche von Stieltjes Stieltjes' the- 
ory is full of important ideas. In Chapter V of [105] Stieltjes gives a discussion on 
the convergence of sums of the form 

/i(*) + / 2 (*) + ••• + /»(*), 

where fi(z) are analytic functions on the open unit disk Cr with center at the 
origin and radius R. He proves a result which was later also proved by Giuseppe 
Vitali in 1903 [118]: 

Theorem (Stieltjes- Vitali). Let f n be a sequence of analytic functions on a non- 
empty connected open set Q of the complex plane. If f n is uniformly bounded on 
compact sets of O and if f n converges on a subset E C O that has an accumulation 
point in Vt, then f n converges uniformly on every compact subset ofVt. 

Paul Montel refers to this theorem as Stieltjes' theorem [67] and others refer to 
it as Vitali's theorem. This result is very convenient in the study of convergence of 
continued fractions because quite often one is dealing with rational fractions and 
one may be able to prove convergence on a set E that is far enough away from 
the poles of the rational fraction. The Stieltjes-Vitali theorem then allows one to 
extend the asymptotic result to hold everywhere except at the set containing all 
the poles. 

The continued fraction (1.2) was later studied by Van Vleck [116] for bk arbitrary 



IMPACT OF STIELTJES' WORK 



5 



then not necessarily supported on [0, oo) and these continued fractions are then 
closely related to Stieltjes integrals over (—00, 00). The complete extension is due 
to Hamburger [37]. Van Vleck [115] and Pringsheim [83] [84] have also given an 
extension to complex coefficients. 

For some good expositions on continued fraction we refer to the books by Perron 
[78], Wall [119], Jones and Thron [46] and Lorentzen and Waadeland [58]. 



2.1 The Stieltjes moment problem. 

In his fundamental work [105, §24] Stieltjes introduced the following problem: 
given an infinite sequence (k = 0, 1, 2, . . . ), find a distribution of mass (a positive 
measure y) on the semi-infinite interval [0, 00) such that 



Of course such a measure will not always exist for any sequence Hk and if such a 
measure exists, then it need not be unique. The Stieltjes moment problem therefore 
has two parts 

(1) find necessary and/or sufficient conditions for the existence of a solution of 
the moment problem on [0, 00), 

(2) find necessary and/or sufficient conditions for the uniqueness of the solution 
of the moment problem on [0, 00). 

Chebyshev had previously investigated integrals and sums of the form 



where w(t) is a positive weight function and Wi are positive weights. Stieltjes 
integrals cover both cases and give a unified approach to the theory. Chebyshev 
did not investigate a moment problem, but was interested when a given sequence 
of moments determines the function w(x) or the weights Wi uniquely. His work 
and the work of his student Markov is very relevant, but Stieltjes apparently was 
unaware of it. See Krein [53] for some history related to the work of Chebyshev and 
Markov. Nevertheless Stieltjes' introduction of the moment problem is still regarded 
as an important mathematical achievement. The reason for the introduction of 
this moment problem is a close connection between ^'-fractions or J-fractions and 
infinite series. If we make a formal expansion of the function 



2. Moment Problems 







which is known as the Stieltjes transform of the measure //, then we find 




This series does not always converge and should be considered as an asymptotic 
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continued fraction of the form (1.1) or (1.2). The n-th approximant of the J- 
fraction has the property that the first 2n terms in the expansion 



J_ Pn-l( X ) 

ai Pn(x) 



oo 

E 

fc=0 



(-1) 



k ™k 
fc+i 



x 



agree with those of the expansion of S(/i; x) i.e., mu = H-k for k = 0, 1, . . . , 2n — 1. 
This rational function is therefore a (diagonal) Pade approximant for S(fi; x). If the 
infinite series is given, then the continued fraction is completely known whenever 
the measure \x is known, provided the continued fraction converges. 

Stieltjes gave necessary and sufficient conditions for the existence of a solution 
of the Stieltjes moment problem: 

Theorem. // the Hankel determinants satisfy 



(2.1) 



and 



tn 



tl2 



tin 
tl2n 



> 0, n£l 



(2.2) 



fJ>l 
^2 



^2 
^3 



l^n+2 



P-n+2 
^2n+l 



> 0, n E N, 



then there exists a solution of the Stieltjes moment problem. 

If the moment problem has a unique solution then the moment problem is de- 
terminate. If there exist at least two solutions then the moment problem is inde- 
terminate. Other terminology is also in use: determined/indetermined and deter- 
mined/undetermined. Any convex combination of two solutions is another solution, 
hence in case of an indeterminate moment problem there will always be an infinite 
number of solutions. Stieltjes gave explicit examples of indeterminate moment prob- 
lems (see also Section 5.3) and he showed that a moment problem is determinate if 
and only if the corresponding continued fraction (1.1) converges for every z in the 
complex plane, except for z real and negative. A necessary and sufficient condition 
for a determinate moment problem is the divergence of the series Yln°=i c n where 
c n are the coefficients of the S- fraction (1.1). In case of an indeterminate moment 
problem Stieltjes constructs two solutions as follows: let P n (z)/Q n (z) be the n-th 
convergent of the continued fraction (1.1), then the limits 

lim P 2n {z) =p{z), lim P 2n+ i(z) =pi{z), 

n — >oo n — >oo 

lim Q 2n {z) = q(z), lim Q 2n +i(z) = qi(z), 



exist, where p,pi,q,qi are entire functions satisfying 
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Stieltjes then shows that 



p( z ) _ r k Pi(z) _ s | s k 



The poles x k , yu (k = 1, 2, ... ) are all real and positive and the residues r k , s k are 
all positive: this follows because the zeros of the numerator polynomials interlace 
with the zeros of the numerator polynomials and because all these zeros are real 
and negative. These limits can thus be expressed as a Stieltjes integral 

p(z) f°° dn{t) Pl (z) r°° dm(t) 



q(z) J z + t qi(z) J z + t 

and both \i and \i\ are solutions of the moment problem with remarkable extremal 
properties. This is one instance where it is clear why Stieltjes introduced the 
concept of a Stieltjes integral. 

Not much work on the Stieltjes moment problem was done after Stieltjes' death. 
One exception is G. H. Hardy [38] who considered the moments of a weight function 
w(x) on [0, oo) with restricted behaviour at infinity: 

r°° 

I w(x)e k ^ dx < oo, 
Jo 

for a positive value of k. He shows that the Stieltjes moment problem is then always 
determinate and constructs the density from the series 

(2n)\ ' 

Hardy's proof avoids the use of continued fractions. 

2.2 Other moment problems. 

Nothing new happened until 1920 when Hamburger [37] extended Stieltjes' mo- 
ment problem by allowing the solution to be a measure on the whole real line 
instead of the positive interval [0, oo). The extension seems straightforward but 
the analysis is more complicated because the coefficients of the continued fraction 
(1.1) may become negative or vanish. Hamburger showed, using continued fraction 
techniques, that a necessary and sufficient condition for the existence of a solution 
of the Hamburger moment problem is the positivity of the Hankel determinants 
(2.1). He also shows that a Hamburger moment problem may be indeterminate 
while the Stieltjes moment problem with the same moments is determinate. 

Nevanlinna [70] introduced techniques of modern function theory to investigate 
moment problems without using continued fractions. His work is important because 
of the notion of extremal solutions, which were first studied by him. M. Riesz [85] 
[86] gave a close connection between the density of polynomials in L 2 -spaces and 
moment problems: 

Theorem. Let fx be a positive measure on (—00,00). If the Hamburger moment 
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are dense in L 2 (/i). If the Hamburger moment problem is indeterminate then the 
polynomials are dense in L 2 (ij) if and only if ji is a Nevanlinna extremal measure. 

Berg and Thill [9] have recently pointed out that this connection is not any 
longer valid in higher dimensions by showing that there exist rotation invariant 
measures // on 1^,(1 > 1 for which the moment problem is determinate but for 
which polynomials are not dense in L 2 (fx). 

In 1923 Hausdorff [39] studied the moment problem for measures on a finite 
interval [a, b] . The Hausdorff moment problem is always determinate and condi- 
tions for the existence of a solution can be given in terms of completely monotonic 
sequences. The moment problem is closely related to quadratic forms of infinitely 
many variables and operators in Hilbert space, as became clear from the work of 
Carleman [12] [13] and Stone [107]. Carleman established the following sufficient 
condition for a determinate moment problem: 



-l/2k 
\i 2k ' =00. 



k=l 

This is still the most general sufficient condition. Karlin and his collaborators [50] 
[51] have approached the moment problem through the geometry of convex sets and 
have shown that many results can be interpreted in this geometrical setting. Let 
me mention here that one can find excellent treatments of the moment problem in 
the monograph of Shohat and Tamarkin [92] and the book of Akhiezer [3] . Also of 
interest is the monograph by Krein and Nudelman [54]. 

2.3 Recent extensions of the moment problem. 

The most recent extension of the moment problem is to consider a doubly infinite 
sequence fi n (n G Z) and to find a positive measure fx on (— oo, oo) such that 



/oo 
x n d/j,(x), n G Z. 
-oo 



Such a moment problem is known as a strong moment problem. The strong Stieltjes 
moment problem was posed and solved by Jones, Thron and Waadeland in 1980 
[48] and again the solution is given in terms of the positivity of certain Hankel 
determinants. These authors again use continued fractions, but instead of the S- 
and J-fractions encountered by Stieltjes and Hamburger, one deals with another 
kind of fraction known as a T-fraction. The strong Hamburger moment problem 
was handled by Jones, Thron and Njastad in 1984 [47]. Njastad [71] gave another 

(k) 

extension, known as the extended moment problem: given p sequences /i„ [n = 
1, 2, . . . ; 1 < k < p) and p real numbers a±, . . . , a p , does there exist a positive 
measure \i on the real line such that 

The solution is again given in terms of positive definiteness of a certain functional. 
Orthogonal polynomials play an important role in the Stieltjes and Hamburger 
moment problem; for the strong moment problem a similar important role is played 
by orthogonal Laurent polynomials and for the extended moment problem one 
deals with orthogonal rational functions. The first place where orthogonal Laurent 
polynomials are considered seems to be a paper by Pastro [75], where an explicit 
example of the orthogonal Laurent polynomials with respect to the Stieltjes- Wigert 
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3. Electrostatic Interpretation of Zeros 

Stieltjes gave a very interesting interpretation of the zeros of Jacobi, Laguerre and 
Hermite polynomials in terms of a problem of electrostatic equilibrium. Suppose 
n unit charges at points 

•£1 5 >^2? • • • ? are distributed in the (possibly infinite) 
interval (a, b). The expression 

D(X!,X 2 , ■ ■ ■ ,X n ) = ] j \Xi-Xj\ 
l<i<j'<n 

is known as the discriminant of x±,X2,--- ,x n . If the charges repell each other 
according to the law of logaritmic potential, then 



-\ogD(x 1 ,x 2 ,... ,x n ) = log- 



l<i<j<n 



is the energy of the system of electrostatic charges and the minimum of this ex- 
pression gives the electrostatic equilibrium. The points x±,X2,... ,x n where the 
minimum is obtained are the places where the charges will settle down. Stielt- 
jes observed that these points are closely related to zeros of classical orthogonal 
polynomials. 

3.1 Jacobi polynomials. 

Suppose the n unit charges are distributed in [—1, 1] and that we add two extra 
charges at the endpoints, a charge p > at +1 and a charge q > at — 1. Each 
of the unit charges interacts with the charges at ±1 and therefore the electrostatic 
energy becomes 

n ^ n \ 

(3.1) L = - \ogD n (x 1 ,x 2 , ■ ,x n ) +p^log— — — j- + q ^ log — — - — -. 

i=l i=l 

Stieltjes then proved the following result [97] [98] [100] 

Theorem. The expression (3.1) becomes a minimum when x±, X2, • • • , x n are the 
zeros of the Jacobi polynomial Pn V l,2q 

Proof:. It is clear that for the minimum all the X{ are distinct and different from 
±1. For a minimum we need dL/dxk = (1 < k < n) so that we have the system 
of equations 

£— ^ ^T-^T = ' l<k<n. 

X{ Xfc Xk 1 Xk -\- 1 

1=1 

If we introduce the polynomial 

p n (x) = (X - Xij(x -X 2 )---{X- X n ), 

then this is equivalent with 



10 



WALTER VAN ASSCHE 



This means that the polynomial 

(1 - x 2 )p'^x) + 2[q-p-(p + q)x\p' n {x) 

vanishes at the points Xk and since this polynomial is of degree n it must be a 
multiple of p n (x). The factor is easily obtained by equating the coefficient of x n 
and we have 

(1 - x 2 )p'^{x) + 2[q-p-(p + q)x]p' n (x) = -n[n + 2{p + q) - l]p n {x), 

which is the differential equation for the Jacobi polynomial P^ v 1,2q 1 ' ) (x)/c n , 
where c n is the leading coefficient of the Jacobi polynomial. □ 

Stieltjes also found the minimum value. Hilbert [42] also computed the minimum 
value and Schur [89] treated the case p = q = in detail. Schur's paper then led 
Fekete [26] to define the transfinite diameter of a compact set K (with infinitely 
many points) in the complex plane. Take n points Zi G K (i = 1, 2, . . . , n), and 
put 

d n = maxD(zi, . . . ^n) 1 ^ 2 ), 

ZiEK 

then d n is a decreasing and positive sequence [26] [111, Thm. III. 21 on p. 71]. The 
limit of this sequence is the transfinite diameter of K and is an important quantity 
in logarithmic potential theory (see §3.4). The transfinite diameter thus comes 
directly from Stieltjes' work. 
Consider the function 

' n \ 

JJ|l-x i | a, - 1 |H-x i | 1 '- 1 D n (x u ... ,x n ), 

\i=l J 

then Stieltjes' electrostatic interpretation gives the ^-norm of this function. 
The L^_ 1 -,^^-norm of this function is also very famous and is known as Selberg's 
beta integral [90]. Actually Selberg evaluated a multiple integral over [0, l] n : 

jf 1 ... jT 1 £>„(*!,... , t n ) 2z (j[ tr\i - uy-^j dh... dt n 



n 



T(x + (j - l)z)T(y + (j - l)z)T(jz + 1) 



T(x + y + (n + j-2)z)T(z + l) 

but this integral can easily be transformed to an integral over [—1, l] n which by 
an appropriate choice of the parameters z,x,y becomes the desired L^_ 1 ^-norm. 
This multiple integral has many important applications e.g., in the statistical theory 
of high energy levels (Mehta [64]) but also in the algebraic theory of root systems 
(Macdonald [60]). Aomoto [6] gave an elementary evaluation of Selberg's integral 
and Gustafson [35] computed some g-extensions. Selberg's work was not inspired 
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3.2 Laguerre and Hermite polynomials. 

A similar interpretation exists for the zeros of Laguerre and Hermite polynomials. 
Suppose the n unit charges are distributed in [0, oo) and that we add one extra 
charge p > at the origin. In order to prevent the charges from moving to oo we 
add the extra condition that the centroid satisfies 

1 n 

(3.2) -$>,<K, 

fc=i 

with K a positive number. The energy now is given by the expression 

1 

(3.3) L = -logD„(xi, . . . ,x n ) + pJ^log — . 

fc=i Xk 

Theorem. The expression (3.3) together with the constraint (3.2) has a minimum 
when xi, X2, ■ ■ ■ ,x n are the zeros of the Laguerre polynomial Ln P (c n x), where 
c n = (n + 2p-l)/K. 

If the n unit charges are on (—00,00) and if the moment of inertia satisfies 

1 n 

(3.4) n^ xl ~ L > 

n k=i 

then 

Theorem. The expression — \ogD n (xi,X2,--- ,x n ) with constraint (3.4) becomes 
minimal when x±,X2,.- - ,x n are the zeros of the Hermite polynomial H n (d n x), 
where d n = a/ \n — 1)/2L. 

The proof of both statements is similar to the proof for the Jacobi case, except 
that now we use a Lagrange multiplier to find the constrained minimum. Mehta's 
book on Random Matrices [64] gives an alternative way to prove the results for 
Laguerre and Hermite polynomials. 

In 1945 Siegel [93] reproved the theorem for Laguerre polynomials and applied 
it to improve the arithmetic-geometric mean inequality and to find better bounds 
on algebraic integers. Siegels seems not to have been aware of Stieltjes' work, but 
started from Schur's work [89]. 

3.3 Extensions. 

In [99] Stieltjes generalizes this idea to polynomial solutions of the differential 
equation 

(3.5) A(x)y" + 2B(x)y' + C(x)y = 0, 

where A, B and C are polynomials of degree respectively p+l,p and p — 1. Such a 
differential equation is known as a Lame equation in algebraic form. In 1878 Heine 
asserted that when A and B are given, there are in general exactly ( n+p ~ 1 ) poly- 
nomials C such that the differential equation has a solution which is a polynomial 
of degree n. Stieltjes assumes that 

B(x) = >A r k 

A ( T \ T — ftp ' 
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with rk > and G K. One can then put charges at the points and n unit 
charges at n points x\,X2, ■ ■ ■ , x n on the real line. Stieltjes then shows that there 
are exactly ( n+ ^ _1 ) positions of electrostatic equilibrium, each corresponding to 
one particular distribution of the n charges in the p intervals [a,k, ajfc+i] (0 < k < p), 
and these charges are then at the points x\, x%, ■ ■ ■ ,x n which are the n zeros of 
the polynomial solution of the differential equation. This result is now known as 
the Heine-Stieltjes theorem [110, Theorem 6.8 on p. 151]. The conditions imposed 
by Stieltjes have been weakened by Van Vleck [114] and Polya [79]. Polya allowed 
the zeros of A to be complex and showed that the zeros of the polynomial solution 
of the differential equation will all belong to the convex hull of {ao, . . . , a p }. The 
location of the zeros of the polynomial solution is still under investigation now 
and interesting results and applications to certain problems in physics and fluid 
mechanics are discussed in [4] [5] [123]. 

Recently Forrester and Rogers [27] and Hendriksen and van Rossum [40] have 
allowed the n unit charges to move into the complex plane. Forrester and Rogers 
consider a system of 2n particles of unit charge confined to a circle in the complex 
plane, say at the points e ldj and e~ l9j (1 < j < n). At 9 = (i.e., at the point 
z = 1) a particle of charge q is fixed and ai 9 = n (z = — 1) a particle of charge p. 
The energy of the system is now given by 

In 2n 

(3.6) L = —q ^ log 1 1 — e i6k \ — p ^ log 1 1 + e i9k \ — ^ log |e lSfc - e^' |, 

fc=l k=l l<k<j<2n 

where 

(3.7) < 6j < 7v, 6 j + 9 n+j =2n, 1 < j < n. 

Theorem (Forrester and Rogers). The minimum of L given in (3.6) subject 
to the constraints (3.7) occurs when 9j are the zeros of the trigonometric Jacobi 

polynomial Pn 2 ' P 2 \cos9). 

Forrester and Rogers also consider crystal lattice structures in which n2 m par- 
ticles of unit charge and 2 m particles of charge q are distributed on the unit circle, 
with one of the q charges fixed at 9 = 0. If one requires that between every two q 
charges there are n unit charges then the equilibrium position of the n2 m particles 

of unit charge occurs at the zeros of the Jacobi polynomial P^l ' 2 (cos 2 m 9) when 

n is even and at the zeros of P( n liy 2 ( cos 2 m 6>) when n is odd. The equilibrium 
position of the 2 m — 1 particles of charge q occurs at 9k = ?0 (1 < k < 2 m ). 

Hendriksen and van Rossum [40] have considered situations where other special 
polynomials come into play. Suppose a > and that there is a charge (a + l)/2 
at the origin and a charge (c — a)/2 at the point 1/a. If a — > oo one obtains a 
generalized dipole at the origin. Suppose now that there are n unit charges at 
points zi,Z2,--. ,z n in the complex plane, then the electrostatic equilibrium in 
this generalized dipole field is obtained when z\ , . . . , z n are the zeros of the Bessel 
polynomial 2-Fo( — n, c + n;x). Similar results can be obtained on so-called m-stars 
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Suppose that positive charges q are placed at the endpoints p = r of S m and a 
charge p > is placed at the origin. If the points zi, . . . ,z n (n > m) in the 
complex plane all have a unit charge, then the electrostatic equilibrium (assuming 
rotational symmetry) is obtained by choosing z\ , . . . ,z n to be the zeros of the 
polynomial f n of degree n that is a solution of the differential equation 

(r m - z m )zy" - 2[(p + qm)z m - pr m ]y' = -n(n - 1 + 2p + 2qm)z m - 1 y. 

For particular choices of the parameters p, r, m one then obtains well known (or- 
thogonal) polynomials. 

3.4 Logarithmic potential theory. 

Suppose that we normalize the electrostatic problem on [—1, 1] in such a way 
that the total charge is equal to 1. The n charges then are equal to l/(n + p + q) 
and the charges at 1 and — 1 become respectively p/(n + p + q) and q/(n + p + q). 
What happens if the number of particles n increases? Clearly the charges at the 
endpoints ±1 become negligible compared to the total charge of the particles inside 
[— 1, 1]. This is the only place where p and q affect the distribution of the zeros, 
therefore it follows that the asymptotic distribution of the charges in (—1, 1) i.e., 

the asymptotic distribution of the zeros of Jacobi polynomials Pn P ~ 1,2q ~ 1 \x), is 
independent of p and q. By takings = q = 1/4 we deal with Chebyshev polynomials 
of the first kind T n {x) with zeros {cos ( - 2j ~ n 1 - >7r : 1 < j < n}. Let N n (a,b) be the 
number of zeros of T n (x) in [a, b], then 



N n (a,b) 1 
n ' n 

a<cos 2n <b 



/ ldt + o(l) 

J a<cos t-K<b 

1 f b dx 



Therefore the asymptotic distribution of the zeros of Jacobi polynomials is given 
by the arcsin distribution and the relative number of zeros in [a, b] is 

1 f b dx 1 

— / = — ( arcsin b — arcsin a). 

The surprising thing is that this is valid not only for Jacobi polynomials but for 
a very large class of orthogonal polynomials on [—1, 1]. The arcsin distribution is 
actually an extremal measure in logarithmic potential theory. Widom [120] [121] 
and Ullman [112] were probably the first to connect logarithmic potential theory and 
general orthogonal polynomials, even though some aspects such as the transfinite 
diameter and conformal mappings had already appeared in earlier work by Szego 
[110, Chapter XVI]. Let K be a compact set in C and denote by Qk be the set of 
all probability measures on K. Define for \x e Ok the logarithmic energy by 



!(»)= f f log-: — 1 — rdu(x)du(y), 
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then there exists a unique measure \ik £ ^k such that 

I(jm k ) = min 

and this measure is the equilibrium measure (see e.g. [111]). When if = [—1,1] 
then the equilibrium measure turns out to be the arcsin measure. The capacity of 
the compact set if is given by 

cap(if) = e- I{ ^ K \ 
and the capacity of a Borel set B e B is defined as 

cap(£?) = sup cap (if), 

K<ZB, K compact 

(the capacity of B is allowed to be oo). Szego [108] showed that the capacity of a 
compact set if is the same as the transfinite diameter of this set, which we defined 
earlier. The following result concerning the asymptotic distribution of zeros of 
orthogonal polynomials is known (see e.g. [94]): 

Theorem. Let fx be a probability measure on a compact set if C K such that 

inf cap(B) = cap(K), 

M (B) = 1,B6B 

where B are the Borel subsets of if , and suppose that Xk, n (1 < k <n) are the zeros 
of the orthogonal polynomial of degree n for the measure fx. Then 



1 n f 
lim - V] f(x k , n ) = / f{t)dfi K 

n^oo n ^ J K 



(t) 



holds for every continuous function f on K . 

When K = [—1,1] then the conditions hold when fi is absolutely continuous 
on (—1,1) with filix) > almost everywhere (in Lebesgue sense). This includes 
all Jacobi weights. A very detailed account of logarithmic potential theory and 
orthogonal polynomials can be found in a forthcoming book by H. Stahl and V. 
Totik [94]. 

There is a similar generalization of the electrostatic interpretation of the zeros 
of Laguerre and Hermite polynomials. This time we need to introduce the energy 
of a measure in an external field /. If if is a closed set in the complex plane C and 
if the field / : K — > [0, oo) is admissible i.e., 

(1) / is upper semi-continuous, 

(2) the set {z G K : f(z) > 0} has positive capacity (oo is allowed), 

(3) if K is unbounded then zf(z) — > as \z\ — > oo (z E if), 

then we define the energy integral in the field / as 

I f(l J ) = [ [ l °Sr~ — \dfi(x)dfi(y)-2 I log f(x) dfi(x). 

JkJk \ x ~y\ JK 

Again there exists a unique measure fx / such that 

If(fUf) = min I f (fi), 

and this measure is the equilibrium measure in the external field / [34] [65]. The 
following result generalizes the electrostatic interpretation of the zeros of Hermite 
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Theorem. Suppose that Xf~, n (1 < k < n) are the zeros of the n-th degree orthog- 
onal polynomial with weight function w(x) on (—00, 00). Suppose that there exists 
a positive and increasing sequence c n such that 

(3.8) lim wM 1/n = /(i), 16I, 

uniformly on every closed interval, with f an admissible field, then 



n^oo fl 

k=l 

for every bounded and continuous function g. 

Again the asymptotic distribution of the (contracted) zeros of orthogonal poly- 
nomials does not depend on the exact magnitude of the weight function w, but only 
on the asymptotic behaviour given in (3.8). When w(x) = e~l x l - the so-called 
Freud weights — then one can take c n = c(a)n 1 ^ a with 



c(a) 



r(f) 



to find that f{x) = e l c ( a ) x l . The corresponding equilibrium measure \if has 
support on [—1, 1] and is absolutely continuous with weight function 

1 f 1 dy a 

7T J\ t \ y/y 2 - t 2 

This is now known as the Nevai-Ullman weight. Notice that the logarithm of the 
external field is the mathematical counterpart of the constraints (3.2) and (3.4) for 
the Laguerre and Hermite polynomials. 

The fascinating aspects of logaritmic potential theory and zeros of orthogonal 
polynomials are very much inspired by Stieltjes' observation that the zeros of Ja- 
cobi, Laguerre and Hermite polynomials actually solve an equilibrium problem in 
electrostatics. 



4. Markov-Stieltjes Inequalities 

In his paper [95] Stieltjes generalized the Gaussian quadrature formula, which 
Gauss gave for the zeros of Legendre polynomials, to general weight functions on 
an interval [a, b]. E. B. Christoffel had given this generalization already seven years 
earlier [19] [30], but Stieltjes' paper is the first that makes a study of the convergence 
of the quadrature formula. The Gaussian quadrature formula approximates the 
integral 

(4.1) / tt(x) d/j,(x), 

J a 

by appropriately summing n function evaluations 



(4.2) 



n 
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This formula has maximal accuracy 2n — 1 i.e., the sum is equal to the integral for 
all polynomials of degree at most 2n — 1, when the quadrature nodes are the zeros 
Xj,n (1 < J ' < n) of the orthogonal polynomial p n (x) of degree n with orthogonality 
measure //, and the quadrature weights A J>n (1 < j < n) are given by 

A . _ - 1 1 



Q'n+lP n {%j,n)Pn+l {%j,n) ^nP n {%j,n)Pn—l {%j,n) 

where we have used the recurrence relation (1.3). These weights are known as the 
Christoffel numbers and have important properties. One of the most important 
properties is their positivity, which follows easily from 



X 3,n = \ ^Pl{Xj,n) 
lfc=0 

Stieltjes gives another remarkable property, namely 

(4.3) ^ Xj, n < / djj(x) = /j,[a,x k , n ) < n[a,x k , n ] < J^ A J>- 

3 = 1 Ja 3 = 1 

Stieltjes was unaware that Chebyshev had already conjectured these inequalities in 
[15] and that Chebyshev's student A. A. Markov had proved them in [62]. Markov's 
paper appeared in 1884, the same year as Stieltjes' paper [95], but in [96] Stieltjes 
kindly acknowledges Markov to be the first to prove the inequalities. He says 
however that his proof is independent of Markov's proof since [95] was submitted 
in May 1884 whereas Markov's paper arrived at the library in September 1884. 
Szego [110] gives three proofs of (4.3), combining the proofs of Stieltjes and Markov. 
The inequalities (4.3) are nowadays known as the Markov-Stieltjes inequalities. A 
related set of inequalities was proved by K. Posse [80] [81]. If / : (a, b) — > 1R. is such 
that (x) > for every x G (a, Xk, n ] (j = 0, 1, 2, . . . , 2n — 1) then 

J2Xj,nf(Xj, n ) < / f(x)dfl(x) < J2Xj,nf(Xj, n ). 

3=1 Ja 3=1 

Stieltjes [96] [106] also gives other inequalities for the Christoffel numbers e.g., 

k-1 k k 

3=1 3=1 3=1 

Stieltjes used the Markov-Stieltjes inequalities to show that the sum (4.2) takes 
the form of a Riemann-Stieltjes sum for the integral (4.1), which makes Stieltjes 
the first to prove convergence of Gaussian quadrature for continuous functions. If 
zeC\[fl,i] then the function f(x) = l/(z — x) is continuous on [a, b] and hence 
the Gaussian quadrature applied to / converges to the Stieltjes transform of the 
orthogonality measure \x. The Gaussian quadrature formula for this function / 
is a rational function of z and is exactly the n-th approximant (n-th convergent) 
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quadrature formula for / therefore gives an important result of Markov regarding 
the convergence of the diagonal in the Pade table for the Stieltjes transform of a 
positive measure. 

Another important application of the Markov-Stieltjes inequalities is a necessary 
and sufficient condition for determinacy of the moment problem: if 



oo 

n=0 

for every real x which is not a point of discontinuity of fi, then the moment prob- 
lem for \i is determinate. The Markov-Stieltjes inequalities are also very useful 
for estimations of the rate of convergence of the Gaussian quadrature formula; the 
Posse-Markov-Stieltjes inequalities even give results for singular integrands (Lubin- 
sky and Rabinowitz [59]). The estimation of the distance between two succesive 
zeros of orthogonal (and also quasi-orthogonal) polynomials can also be done using 
these inequalities. From (4.3) one finds 

(4.5) X jin < fi(xj + i, n ,Xj- ltn ), 

which allowed Stieltjes to deduce that A J;n tends to zero when n — > oo whenever 
the behaviour of Xj + i tn — Xj-i jTl is known in terms of the measure \i. Stieltjes 
gave the result for Legendre polynomials. Nevai [68, p. 21] used the bound (4.5) 
to show that for measures ji with compact support such that for every e > the 
set supp(/u) fl (x — e, x + e) is an infinite set, there exists a sequence of integers k n 
(n = 1, 2, . . . ) such that 

lim x kn , n = x, lim A fc „, n = 0. 

n — >oo n — >oo 

This shows that the zeros are dense in the derived set of supp(/u) and that the 
corresponding Christoffel numbers tend to zero. For the isolated points in supp(ju) 
Nevai [68, p. 156] used (4.3) to show that 

lim limsup X k , n = l*({x}), 

\x-Xk,n\<£ 

for every x G R. Freud [28, p. Ill] shows that for two consecutive zeros in an 
interval [c, d] for which 

0<m <^>lM<M, x,ye\c,d], 
y-x ' 

one has 

ci _ C2 

_ Xj+l, n %j,n _ 7 

n n 
where c\ , C2 are positive constants. This is a slight extension of a result by Erdos 
and Turan [24]. Nevai [68, p. 164] generalizes this result by allowing fi' to have an 
algebraic singularity inside supp(/u). If supp(/u) is compact, A C supp(/f), t G A° 
(the interior of the set A) and if [x is absolutely continuous in A with 

c\\x — ty < /i'(x) < ci\x — t| 7 , 7 > —1, 

then 

Si < - < £i 

_ ^fe.n -^fc— l,n — > 

n n 
whenever x k n G Ai with A 1 a closed subset of A°. The Markov-Stieltjes inequali- 
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5. Special Polynomials 

5.1 Legendre polynomials. 

Stieltjes wrote a number of papers directly related to the Legendre polynomials 
P n (x) for which 



J P n (x)Pm(x) dx = 0, m^n. 



He always uses the notation X n but here we will adopt the notation P n which is 
nowadays standard. In [100] he uses the electrostatic interpretation of the zeros 
of Jacobi polynomials to obtain monotonicity properties of the zeros of Jacobi 
polynomials as a function of the parameters, and from this one easily finds bounds 
for the zeros xi jH > X2, n > • • • > x njn of the Legendre polynomials P n (x) = 

Pn ' ; (x) in terms of the zeros of the Jacobi polynomials P n 2 ' 2 (x) and P n 2 ' 2 (x) 
giving 

2kn (2k - 1W 

cos- — — <x fc , n <cos— — — — , l<k<n. 
2n + 1 2n + 1 

These bounds were already given by Bruns [11] in 1881 and Stieltjes does refer to 
Brans' result, but Stieltjes goes on and shows that by using the zeros of Chebyshev 

polynomials of the first kind T n (x) = P n 2 ' 2 (x) and of the second kind U n (x) = 

-Pn 2 ' 2 (x) one may find better bounds: 

kn (2/e-lW 

cos < xu n < cos , 1 < k < n/2. 

n+1 ' 2n - - i 

A great deal of work has been done to obtain sharp bounds for zeros of orthogonal 
polynomials. The monotonicity of zeros of orthogonal polynomials depending on 
a parameter is often used. Markov [61] [110, Theorem 6.12.1] gave a very nice 
result concerning the dependence of the zeros on a parameter t which appears in 
the weight function w(x) = w(x;t). Two other methods for obtaining bounds for 
zeros of orthogonal polynomials are the Sturm comparison theorem [110, §6.3] [57] 
for solutions of Sturm-Liouville differential equations and the Hellman-Feynman 
theorem [45] of quantum chemistry. See also [44] for results on the monotonicity of 
zeros of orthogonal polynomials. 

Stieltjes made a very important contribution concerning the asymptotic be- 
haviour of Legendre polynomials. In 1878 Darboux [23] gave an asymptotic series 
for the Legendre polynomial: 

"'~ ! 1-3- ••(2/,- 1) 



P n (cos0) =2a n ^ a fc (2w _ 1)(2w _ 3) ^ (2w _ 2fc + 1) 



cos\(n-k+±)6-(k+±)ic/2] 

—± 2> \ + 0(n~ m - * , < 6 < tt, 

(2 sm^) fc+ 2 

which generalizes an asymptotic formula given by Laplace (when m = 1). Here 

1 • 3 • • • (2k - 1) 
«o = l, a k = ^ . 

The problem with this formula is that there is no closed expression or a bound 
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sense when ^ < 9 < but it converges to 2P n (cos6>) rather than P n (cos6>) (this 
"paradox" was first pointed out by Olver [72] ) . This is probably the easiest example 
showing that asymptotic expansions need not converge to the function that they 
approximate. The reason why things go wrong here is that the formula is obtained 
by the so-called method of Darboux which consists of obtaining asymptotic results 
of a sequence by carefully examining the singularities on the circle of convergence 
of the generating function. The generating function of Legendre polynomials has 
two singularities on the circle of convergence, and at each singularity one picks 
up information on P n (cos9). This is probably the reason why the convergence of 
the infinite series is to 2P n (cos6>) rather than P n (cos6>). Stieltjes' generalization 
of Laplace's asymptotic formula for the Legendre polynomials does not suffer from 
either problem. Stieltjes' asymptotic expansion is [102] [103] 



P n (cos9) 



2 n nl 



n 3 • 5 • • • (2n + 1) 



(2 sin (9) ^ 1 3 



where 



i 2 -3 2 ---(2fc-i; 

b = 1, b k = 



2 fc fc!(2n + 3)(2n + 5) • • • (2n + 2k + 1) ' 
and the error R m (9) is bounded by 

4 2 n n! M 



\Rm{0)\<b r , 



7r3-5---(2n+l) (2sin6») m +5 



where 

f l/cos0, if sin 2 # < ± 
M=\ ' o ~ ? 

I 2 sin 9, ifsin 2 #>±. 

This asymptotic expansion converges in the ordinary sense when ^ < 9 < ^ and 
it converges to ^^(cos^). Combined with Mehler's asymptotic formula 

q 

lim P n (cos -) = J (9), 

one then finds an asymptotic series for the Bessel function Jo already obtained by 
Poisson, but now with a bound on the error. Stieltjes also uses the asymptotic series 
to obtain approximations of the zeros of the Legendre polynomials. The asymptotic 
theory of orthogonal polynomials (in particular classical orthogonal polynomials) is 
very well developped nowadays, at least for orthogonality on a finite interval. Szego 
has a very nice chapter on the asymptotic properties of the classical polynomials 
[110, Chapter VIII] and that book is still a very good source for asymptotic formulas 
for Jacobi, Laguerre and Hermite polynomials. 

A third contribution of Stieltjes involving Legendre polynomials is his work on 
Legendre functions of the second kind [104]. The Legendre function of the second 
kind can be defined by 

(5.1) Qn(x) = l [ ^-dy, xeC\[-l,l], 
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so that 

(5-2) Q n (x) = \p n {x) log (|±±) - ^(x), 

where Pj l 1 } 1 (x) is the associated Legendre polynomial of degree n — 1. The integral 
representation cannot be used to define Q n (x) for x G [—1,1] but by taking the 
appropriate limit, or the appropriate branch of the logarithm, one can use (5.2) to 
define Q n (x) for — 1 < x < 1. Hermite [41] had studied the zeros of Q n on [—1, 1] by 
making some changes of variables. Stieltjes works directly with Q n (x) as a function 
of the real variable x and shows that Q n (x) has n + 1 zeros in (—1, 1) which interlace 
with the zeros of the Legendre polynomial P n (x). He also shows that there can be 
no zeros outside [—1, 1] by using a simple property of Stieltjes transforms of positive 
weight functions. Some of these results can easily be generalized to functions of the 
second kind corresponding to general orthogonal polynomials [113]. 

5.2 Stieltjes polynomials. 

In his last letter to Hermite [8, vol. II, pp. 439-441] Stieltjes considers the Le- 
gendre functions of the second kind (5.1) and observes that 

1 „ , x Q>~\ Q-2 

E n+1 (z) + — + + , 



Qn(z) "' Z Z 

where E n +i(z) is a polynomial of degree n + 1. This polynomial is now known as 
the Stieltjes polynomial 1 . Stieltjes gives the remarkable property 



i 

k 



P n (x)E n+1 (x)x k dx = 0, < k < n, 

which essentially means that E n +i(x) is orthogonal to all polynomials of degree 
less than or equal to n with respect to the oscillating weight function P n (x) on 
[—1,1]. One may now wonder which properties of ordinary orthogonal polynomials 
are still valid for E n+ i(x) and Stieltjes conjectures that the zeros of E n+ \{x) are 
real, simple and belong to [—1, 1] and that they interlace with the zeros of P n (x). 
These conjectures were later proved by Szego [109]. Szego also extended the idea 
to ultraspherical weights by considering the functions of the second kind 



2r(/i + |) i_i z-x 

where P£(x) is an ultraspherical polynomial of degree n. One can then find 

1 



q%(z) n+1K ' z z* 
where E£ +1 (z) is a polynomial of degree n + 1. Szego shows that 



l 

2\u- 



(1 - x z y-*P%{x)E^ +1 (x)x k dx = 0, < k < n, 
i 



lr The polynomial solutions of a Lame differential equation (3.5) are also known as Stieltjes 
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thus generalizing the orthogonality of Stieltjes polynomials. The properties of the 
zeros of E^ +1 (x) depend on the value of the parameter \i. If < ji < 2 then the 
zeros of E^ +1 (x) are in [—1, 1], they are real and simple and they interlace with 
the zeros of P£(x). When \x < then some of the zeros are outside [—1,1] and 
Monegato [66] has made some computations showing that for [i > 4.5 there can be 
complex zeros, depending on the degree n. More precise numerical information for 
Gegenbauer weights as well as for Jacobi weights has been obtained by Gautschi 
and Notaris [32]. The construction of Stieltjes and Szego can be generalized by 
considering a positive measure \x on R. Suppose that p n (x; fi) (n = 0, 1, 2, . . . ) are 
the orthogonal polynomials with respect to the measure fi, then the functions of 
the second kind are 

q n (z;/j) = / — — - — d/j,(x), 



z — x 



and these are defined for z G C\supp(/t). Define the (general) Stieltjes polynomial 

E n+1 {z;n) by 

E n+1 (z;n) H 1 = 1 , 



q n {z;H) z z z 

then one always has 

J p n (x; n)E n+ i(x; fi)x k dfi(x) = 0, < k < n. 

These Stieltjes polynomials turn out to have some importance in constructing 
an optimal pair (A, B) of quadrature formulas. Suppose we start with a quadrature 
formula A with n nodes and a quadrature formula B with m nodes (m > n). In 
order to compute the error of formula A one often assumes that the difference of 
the results obtained by using A and B is proportional to the actual error of the 
quadrature formula A. This means that one needs n + m function evaluations to 
compute the error of A. This implies that one has done m extra function evaluations 
which are not used in the evaluation of A itself. Kronrod [55] suggested to extend 
formula B to a formula with n + m nodes in such a way that the accuracy of B is as 
high as possible. For the Legendre weight on [—1, 1] one will find an optimal pair 
(A, B) by taking for A the Gaussian quadrature with nodes equal to the zeros of 
the Legendre polynomial P n (x) and for B a quadrature formula with 2n + 1 nodes 
at the zeros of P n (x) and the zeros of E n+ i(x). The quadrature formula B then 
turns out to give a correct result for all polynomials of degree less than or equal to 
3n + 1 [55, Theorem 6]. 

In 1930 Geronimus [33] slightly changes Stieltjes' idea and considers the Jacobi 
functions of the second kind 

Q (a,(l) {z)= ^ W {1 _ x) a {1+x)f 3 d ^ 

y_i z-x 

with Pi a,f3 \x) the Jacobi polynomial of degree n. Geronimus observes that 

1 C f \ , Cl , C2 
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with S n (z) a polynomial of degree n. Notice the extra factor \J z 1 — 1 in the denom- 
inator on the left hand side. These polynomials satisfy the remarkable property 

J (l-x) a (l+x) f3 P^ f3 \x)S n (x)T k (x)dx = 7 0<k<n, 

and 

J {l-x) a {l + xfP^\x)S n {x)dx = l. 

Here T k (x) is the Chebyshev polynomial of the first kind of degree k. Geronimus 
polynomials can be generalized to other weights on [—1,1]. The interval [—1,1] 
is important because it accounts for the factor \J z 1 — 1 in the definition of the 
Geronimus polynomials. There is a relation between the Geronimus polynomials 
S n (x) and the Stieltjes polynomials E n+ i(x) if one works with a weight function 
on [-1,1]: if 

n+l 

E n+1 (x) = c kj7l T k (x), 

k=0 

(the prime means to divide the first term by two) is the expansion of E n+ i(x) in 
Chebyshev polynomials of the first kind, then 

n 

S n(x) = ^2c k +l, n U k (x) 
k=0 

is the expansion of S n (x) in Chebyshev polynomials of the second kind. 

Stieltjes and Geronimus polynomials and the related Gauss-Kronrod quadrature 
are still being studied and we refer to Gautschi [31], Monegato [66], Peherstorfer 
[76] [77] and Prevost [82] for more information. 

5.3 Stieltjes- Wigert polynomials. 

In his memoir [105, §56] Stieltjes explicitly gives an example of a moment problem 
on [0, 00) which is indeterminate. He shows that 

/ u k u- logu [1 + Xsm(27v log u)] du = v^Fe 1 ^ 
Jo 

is independent of A and therefore the weight functions 

w x (u) = u~ logu [1 + Asin(27rlogw)] , -1 < A < 1 

all have the same moments which implies that this moment problem is indetermi- 
nate. Stieltjes gives the coefficients of the continued fraction (1.1) 

2n + l 

q 2 

(g;i)n 

where q = e -1 / 2 and 
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Both the series J2 c 2n and J2 c 2n+i converge since < q < 1, which agrees with 
the theory worked out by Stieltjes. Later Wigert [122] extended this by considering 
the weight functions 



w k (x) = e fc2lo g 2x 7 < x < oo, 

which for k = 1 reduce to the weight function considered by Stieltjes. If we set 
q = e _1// ( 2fc ) then the orthogonal polynomials are given by 



Pn{x) = j2 (f^v 2 / v +1 

j=0 W ' q)j 



xV. 



and are known as Stieltjes- Wigert polynomials. The moment problem is indeter- 
minate whenever < q < 1, which means that there exist an infinite number of 
measures on [0, oo) with the same moments. Askey [7] indicated that these poly- 
nomials are related to theta functions and shows that the weight function 

£-5/2 

W(x) = - r ; r , < X < OO 

{-x; q)oo {-q/x; q)oo 

has the same moments. This measure arises as a (/-extension of the beta density 
on [0, oo). Chihara [16] [18] has given many more measures which have the same 
moments as the weight function Wk{x) given by Wigert. 

The Stieltjes- Wigert polynomial p n (x) is a (terminating) basic hypergeometric 
series. Such series are of the form J2 c j with Cj + i/cj a rational function of qi for 
a fixed q (for hypergeometric series this ratio is a rational function of j). The first 
set of orthogonal polynomials which are basic hypergeometric series was found by 
Markov in his thesis [63]. Except for a reference in Szego's book [110, §2.9], this work 
was overlooked and seems not to have led to any extensions. Markov's polynomials 
are discrete extensions of Legendre polynomials and basic hypergeometric exten- 
sions of discrete Chebyshev polynomials which are orthogonal on {0, 1,2, .. . , N} 
with respect to the uniform distribution. They are a special case of polynomials 
considered by Hahn, which will be mentioned later. The next basic hypergeomet- 
ric orthogonal polynomials were introduced in 1894, and there were two different 
examples that year. These are the Stieltjes- Wigert polynomials (with q = e -1 / 2 ) 
given by Stieltjes and the continuous g-Hermite polynomials given by Rogers [87]. 
Both are basic hypergeometric extensions of Hermite polynomials but of a com- 
pletely different nature. Those of Rogers are orthogonal on [—1, 1] with respect to 
the weight function 



w[x) -^o [l-2(2x 2 -l)^+^] 



A number of other examples were found before Hahn [36] considered the following 
problem: find all sets of orthogonal polynomials p n (x) (n = 0, 1, 2, . . . ) such that 

p n +i(x) -p n+1 (qx) 
r„lx) = . 71=11.1.2 
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is again a set of orthogonal polynomials. Earlier it had been shown that if p n (x) 
(n = 0, 1, 2, . . . ) are orthogonal and p' n+1 (x) (n = 0, 1, 2, . . . ) are orthogonal, then 
p n (x) (n = 0,1,2,...) are either Jacobi, Laguerre or Hermite polynomials (after 
a possible change of scale). It is easy to see that the Stieltjes-Wigert polynomials 
are in the Hahn class. The continuous (/-Hermite polynomials of Rogers are not 
in the Hahn class, but their analogous difference operator is a divided difference 
operator. Basic hypergeometric series and orthogonal polynomials which are ter- 
minating basic hypergeometric series are described in detail in the book by Gasper 
and Rahman [29]. All of these polynomials arise in the study of quantum groups 
(see Koornwinder [52] and references there). 

5.4 Orthogonal polynomials related to elliptic functions. 

In Chapter XI of his memoir [105] Stieltjes gives some examples of continued 
fractions and the corresponding moment problem. These examples (except for one) 
had already been worked out in one of his previous papers [101]. The continued 
fractions are for the functions 

poo poo 

F 1 (z,k)= cn(w, k)e~ zu du, F 2 (z,k) = dn(u, k)e~ zu du, 
Jo Jo 

pOO rOO 

F 3 (z,k)= / sn(w, k)e~ zu du, F 4 (z,k) = z sn 2 (w, k)e~ zu du, 
Jo Jo 

which are all Laplace transforms of the Jacobian elliptic functions given by 

, 2tt ^ q n-\ (2n- 1)ttu 

cn(«, k) = cos <p = — ^ cos 2R , 

n=l 

, ,v . 2tt ^ q n ~\ . (2n- lWu 

smu, k) = sin = - — > — - sin , 

V ' 1 r kK ^ 1 - q 2n - 1 2K 

n=l 



dn(u, k) = \Jl-k 2 sin 2 y = ^ + ^ - 



oo 

q U7TU 
COS 



with 



u 



Jo Jl 



< 1 + q 2n K ' 
dff 



V 7 !-^ 2 sin 2 6»' 



and 



y/(l ~X 2 )(1 ~k 2 X 2 ) 

The Chudnovsky's [20, p. 197] pointed out that these continued fractions are some 
of the very rare cases where both the function and its continued fraction expansion 
are known explicitly. There are quite a few cases known when the function is given 
in terms of (basic) hypergeometric series and the numerators and denominators of 
the convergents of the continued fraction are classical orthogonal polynomials (in 
Askey's definition). The three-term recurrence relation then gives the coefficients 
of the J-fraction. The functions Fi(z,k) (i = 1,2,3,4) however are not of (basic) 
hypergeometric type and the corresponding orthogonal polynomials are therefore 
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he obtained S'-fractions for Fi and F 2 and J-fractions for F 3 and F 4 . His method 
consists of decomposing a quadratic form with infinitely many variables as a sum 
of squares: 

oo oo 

^2 E a i+3 X i X 3 = c 0(>0 + ao.l^l + a 0,2^2 H f 

i=0 j=0 

+ c 1 (x 1 + a lj2 x 2 + a 1}3 x 3 H ) 2 

+ c 2 (x 2 + a 2 ,3^3 + a 2 ,4^4 H ) 2 H • 

The coefficients of the J-fraction of 

(-l) n a n 



E 

n=0 



z n+l 



are then determined by the coefficients co,c\, . . . and a^+i (i = 0, 1, 2, . . . ). Such 
a decomposition can easily be made when the function 



x n 



n=0 

satisfies an addition formula of the type 

fix + y) = cof(x)f(y) + Cl h{x)h{x) + c 2 f 2 (x)f 2 (y) + ■■■ 

where f m (x) = 0(x m ), as Rogers [88] pointed out when he was reviewing Stieltjes' 
technique. The addition formulas for the Jacobian elliptic functions then readily 
give the desired continued fractions. The orthogonal polynomials that appear are 
defined by the recurrence relations 

C n+1 (x) = xC n (x) - a n C n _i(», D n+1 (x) = xD n (x) - l3 n D n _i{x), 

with 

a 2n = (2n) 2 k 2 , a 2n+1 = (2n+l) 2 , (3 2n = (2n) 2 , (3 2n+1 = (2n + l) 2 k 2 . 

These polynomials have later been studied in detail by Carlitz [14] . The generating 
function for the orthogonal polynomials satisfies a Lame differential equation 

..1(1 1 1 1 , b-n(n+l)x 



2 \x x — 1 x — a) 4x(x — l)(x — a) 

with n = (in general n is an integer). There exist 2n + l values of the parameter 
b for which this Lame equation has algebraic function solutions. Stieltjes approach 
has been generalized to continued fraction expansions for which the generating 
functions of the numerators and denominators of the convergents satisfy a Lame 
differential equation with n ^ 0. This has been done by the Chudnovsky brothers 
[20, pp. 197-201] [21, §13]. Some of these generalizations have interesting applica- 
tions in number theory: the irrationality and bounds on the measure of irrationality 
of some values of complete elliptic integrals of the third kind can be obtained from 



26 



WALTER VAN ASSCHE 



Acknowledgements 

I never quite realized how much work is needed to analyse Stieltjes' work a 
century after his death. I have spent a lot of time in various libraries and received 
a lot of help from the librarians. I would also like to thank various colleagues 
for suggestions, comments and for pointing out omissions and misinterpretations. 
A sincere word of thanks in particular to Marcel de Bruin, Ted Chihara, Walter 
Gautschi, Tom Koornwinder and Doron Lubinsky. Of course nothing would have 
been possible without the help of Gerrit van Dijk: many thanks for having started 
this whole project. 

References 

1. CEuvres Completes de Thomas Jan Stieltjes, vol. I, P. Noordhoff, Groningen, 1914. 

2. (Euvres Completes de Thomas Jan Stieltjes, vol. II, P. Noordhoff, Groningen, 1918. 

3. N. I. Akhiezer, The Classical Moment Problem, Oliver and Boyd, Edinburgh, 1965. 

4. M. Alam, Zeros of Stieltjes and Van Vleck polynomials, Trans. Amer. Math. Soc. 252 (1979), 
197-204. 

5. A. M. Al-Rashed and N. Zaheer, Zeros of Stieltjes and Van Vleck polynomials and applica- 
tions, J. Math. Anal. Appl. 110 (1985), 327-339. 

6. K. Aomoto, Jacobi polynomials associated with Selberg integrals, SIAM J. Math. Anal. 18 
(1987), 545-549. 

7. R. Askey, Orthogonal polynomials and theta functions, Theta Functions, Bowdoin 1987 (part 
2) (L. Ehrenpreis, R. C. Gunning, eds.), Proceedings of Symposia in Pure Mathematics, 
vol. 49, Amer. Math. Soc, Providence, RI, 1989, pp. 299-321. 

8. B. Baillaud and H. Bourget, Correspondance d'Hermite et de Stieltjes, vol. I-II, Gauthier- 
Villars, Paris, 1905. 

9. C. Berg and M. Thill, Rotation invariant moment problems, Acta Math. 167 (1991), 207- 
227. 

10. C. Brezinski, History of Continued Fractions and Fade Approximants, Springer Series in 
Computational Mathematics, vol. 12, Springer- Verlag, Berlin, 1991. 

11. H. Bruns, Zur Theorie der Kugelfunktionen, J. Reine Angew. Math. 90 (1881), 322-328. 

12. T. Carleman, Sur le probleme des moments, C.R. Acad. Sci. Paris 174 (1922), 1680-1682. 

13. , Sur les equations integrates singulieres a noyau reel et symmetrique, Uppsala Uni- 

versitets Arsskrift, 1923, 228 pp. 

14. L. Carlitz, Some orthogonal polynomials related to elliptic functions, Duke Math. J. 27 
(1960), 443-460. 

15. P. L. Chebyshev, Sur les valeurs limites des integrates, J. Math. (2) 19 (1874), 157-160. 

16. T. S. Chihara, A characterization and a class of distribution functions for the Stieltjes- 
Wigert polynomials, Canad. Math. Bull. 13 (1970), 529-532. 

17. , An Introduction to Orthogonal Polynomials, Gordon and Breach, New York, 1978. 

18. , On generalized Stieltjes- Wigert and related orthogonal polynomials, J. Comput. 

Appl. Math. 5 (1979), 291-297. 

19. E. B. Christoffel, Sur une classe particuliere de fonctions entieres et de fonctions continues, 
Ann. Mat. Pura Appl. (2) 8 (1877), 1-10. 

20. D. V. Chudnovsky and G. V. Chudnovsky, Transcendental methods and theta- functions, 
Theta Functions, Bowdoin 1987 (part 2) (L. Ehrenpreis, R. C. Gunning, eds.), Proceedings 
of Symposia in Pure Mathematics, vol. 49, Amer. Math. Soc, Providence, RI, 1989, pp. 167- 
232. 

21. , Computer assisted number theory with applications, Lecture Notes in Mathematics 

1240, Springer- Verlag, Berlin, 1987, pp. 1-68. 

22. E. Cosserat, Notice sur les travaux scientifiques de Thomas- Jean Stieltjes, Ann. Fac. Sci. 
Tolouse 9 (1895), 1-64. 

23. G. Darboux, Memoire sur I 'approximation des fonctions de tres grands nombres, J. Math. 
(3) 4 (1878), 5-56; 377-416. 

24. P. Erdos and P. Turan, On interpolation II, Ann. of Math. 39 (1938), 703-724. 

c\r t m_ _i rif___7__ .__T___^_____ _t _ m _ t_ _ t_ - _ t_ _ rc /~i i~t A _ _ _i n _• n__'_ rtnn /ir\or\ r\r\rr\r\r\r-n 



IMPACT OF STIELTJES' WORK 



27 



26. M. Fekete, Uber die Verteilung der Wurzeln bei gewissen algebraischen Gleichungen mit 
ganzzahligen Koefficienten, Math. Z. 17 (1923), 228-249. 

27. P. J. Forrester and J. B. Rogers, Electrostatics and the zeros of the classical polynomials, 
SIAM J. Math. Anal. 17 (1986), 461-468. 

28. G. Freud, Orthogonal Polynomials, Akademiai Kiado and Pergamon Press, Budapest and 
Oxford, 1971. 

29. G. Gasper and M. Rahman, Basic Hypergeometric Series, Encyclopedia of Mathematics and 
its Applications, vol. 35, Cambridge University Press, Cambridge, 1990. 

30. W. Gautschi, A survey of Gauss-Christoffel quadrature formulae, E.B. Christoffel: the In- 
fluence of his Work on Mathematics and the Physical Sciences (P. L. Butzer, F. Feher, eds.), 
Birkhauser, Basel, 1981, pp. 72-147. 

31. , Gauss-Kronrod quadrature — a survey, Numerical Methods and Approximation 

Theory III (G. V. Milovanovic, ed.), Faculty of Electronic Engineering, University of Nis, 
1988, pp. 39-66. 

32. W. Gautschi and S. E. Notaris, An algebraic study of Gauss-Kronrod quadrature formulae 
for Jacobi weight functions, Math. Comp. 51 (1988), 231-248. 

33. Ya. L. Geronimus, On a set of polynomials, Ann. of Math. 31 (1930), 681-686. 

34. A. A. Gonchar and E. A. Rakhmanov, Equilibrium measure and the distribution of zeros of 
extremal polynomials, Math. USSR Sb. 53 (1986), 119-130. 

35. R. A. Gustafson, A generalization of Selberg's beta integral, Bull. Amer. Math. Soc. 22 
(1990), 97-105. 

36. W. Hahn, Uber Orthogonalpolynome die q-Differenzengleichungen geniigen, Math. Nachr. 2 
(1949), 4-34. 

37. H. Hamburger, Uber eine Erweiterung des Stieltjesschen Momentenproblems, Math. Ann. 
81 (1920), 235-318; 82 (1921), 120-164; 168-187. 

38. G. H. Hardy, Notes on some points in the integral calculus, XLVI On Stieltjes 'probleme 
des moments', Messenger of Math. 46 (1917), 175-182; XLVII On Stieltjes 'probleme des 
moments' (cont.), Messenger of Math. 47 (1918), 81-88. 

39. F. Hausdorff, Momentprobleme fur ein endliches Intervall, Math. Z. 16 (1923), 220-248. 

40. E. Hendriksen and H. van Rossum, Electrostatic interpretation of zeros, Orthogonal Polyno- 
mials and their Applications (M. Alfaro et al., eds.), Lecture Notes in Mathematics, vol. 1329, 
Springer- Verlag, Berlin, 1988, pp. 241-250. 

41. Ch. Hermite, Sur les racines de la fonction spherique de seconde espe.ce, Ann. Fac. Sci. 
Toulouse 4 (1890), 11-10. 

42. D. Hilbert, Uber die Discriminante der im Endlichen abbrechenden hypergeometrische Rei- 
hen, J. Reine Angew. Math. 103 (1888), 337-345. 

43. , Grundziige einer allgemeinen Theorie der linearen Integralgleichungen, Leipzig, 

1912; Chelsea, New York, 1953. 

44. M. E. H. Ismail and M. E. Muldoon, A discrete approach to monotonicity of zeros of or- 
thogonal polynomials, Trans. Amer. Math. Soc. 323 (1991), 65-78. 

45. M. E. H. Ismail and R. Zhang, On the Hellman-Feynman theorem and the variation of zeros 
of certain special functions, Adv. Appl. Math. 9 (1988), 439-446. 

46. W. B. Jones and W. J. Thron, Continued Fractions: Analytic Theory and Applications, 
Encyclopedia of Mathematics and its applications, vol. 11, Addison- Wesley, Reading, Mas- 
sachusetts, 1980. 

47. W. B. Jones, W. J. Thron and O. Njastad, Orthogonal Laurent polynomials and the strong 
Hamburger moment problem, J. Math. Anal. Appl. 98 (1984), 528-554. 

48. W. B. Jones, W. J. Thron and H. Waadeland, A strong Stieltjes moment problem, Trans. 
Amer. Math. Soc. 261 (1980), 503-528. 

49. S. Karlin and J. McGregor, The differential equations of birth and death processes and the 
Stieltjes moment problem, Trans. Amer. Math. Soc. 85 (1957), 489-546. 

50. S. Karlin and L. S. Shapley, Geometry of Moment Spaces, Memoirs Amer. Math. Soc, 
vol. 12, Providence, R.I., 1953. 

51. S. Karlin and W. J. Studden, TchebychefJ Systems: with Applications in Analysis and Sta- 
tistics, Interscience, New York, 1966. 

52. T. H. Koornwinder, Orthogonal polynomials in connection with quantum groups, Orthogonal 



2<S 



WALTER VAN ASSCHE 



53. M. G. Krein, The ideas of P. L. Chebysheff and A. A. Markov in the theory of limiting 
values of integrals and their further development, Amer. Math. Soc. Transl. Series 2, 12, 
Amer. Math. Soc, Providence, RI, 1959, pp. 1-122. 

54. M. G. Krein and A. A. Nudelman, The Markov Moment Problem and Extremal Properties, 
Translations of mathematical monographs, vol. 50, Amer. Math. Soc, Providence, RI, 1977. 

55. A. S. Kronrod, Nodes and Weights for Quadrature Formulas, Nauka, Moscow, 1964 (Rus- 
sian); Consultants Bureau, New York, 1965. 

56. W. Ledermann and G. E. H. Reuter, Spectral theory for the differential equations of simple 
birth and death processes, Philos. Trans. Roy. Soc. London Ser. A 246 (1954), 321-369. 

57. L. Lorch, Elementary comparison techniques for certain classes of Sturm- Liouville equations, 
Differential equations, Sympos. Univ. Upsaliensis Ann. Quigentesimum Celebrantis, vol. 7, 
Almqvist and Wiksell, Stockholm, 1977, pp. 125-133. 

58. L. Lorentzen and H. Waadeland, Continued Fractions and Some of its Applications, North- 
Holland Studies in Computational Mathematics, North-Holland, 1991 (to appear). 

59. D. S. Lubinsky and P. Rabinowitz, Rates of convergence of Gaussian quadrature for singular 
integrands, Math. Comp. 43 (1984), 219-242. 

60. I. G. Macdonald, Orthogonal polynomials associated with root systems, Orthogonal Polyno- 
mials: Theory and Practice [69], pp. 311-318. 

61. A. A. Markov, Sur les racines de certaines equations, Math. Ann. 27 (1886), 177-182. 

62. , Demonstration de certaines inegalites de M. Tchebycheff, Math. Ann. 24 (1884), 

172-180. 

63. , On Some Applications of Algebraic Continued Fractions, Thesis (in Russian), St. 

Petersburg, 1884, 133 pp.. 

64. M. L. Mehta, Random Matrices and the Statistical Theory of Energy Levels, Academic Press, 
New York, 1967 (2nd edition); 1990 (revised and enlarged 2nd edition). 

65. H. N. Mhaskar and E. B. Saff, Where does the LP -norm of a weighted polynomial live?, 
Trans. Amer. Math. Soc. 303 (1987), 109-124. 

66. G. Monegato, Stieltjes polynomials and related quadrature rules, SIAM Review 24 (1982), 
137-157. 

67. P. Montel, Legons Sur les Families Normales de Fonctions Analytiques et leurs Applications, 
Gauthier-Villars, Paris, 1927. 

68. P. G. Nevai, Orthogonal Polynomials, Memoirs Amer. Math. Soc, vol. 213, Providence, R.I., 
1979. 

69. P. Nevai (ed.), Orthogonal Polynomials: Theory and Practice, NATO-ASI series C, vol. 294, 
Kluwer, Dordrecht, 1990. 

70. R. Nevanlinna, Asymptotische Entwickelungen beschrdnkter Funktionen und das Stieltjess- 
che Momentenproblem, Ann. Acad. Sci. Fenn. Ser. A 18 (1922), no. 5 (52 pp.). 

71. O. Njastad, An extended Hamburger moment problem, Proc Edinburgh Math. Soc. 28 
(1985), 167-183. 

72. F. W. J. Olver, A paradox in asymptotics, SIAM J. Math. Anal. 1 (1970), 533-534. 

73. H. Pade, Sur la representation approchee d'une fonction par des fractions rationelles, Ann. 
Sci. Ecole Norm. Sup. (3) (supplement) 9, 1-93. 

74. , Sur la fraction continue de Stieltjes, C.R. Acad. Sci. Paris 132 (1901), 911-912. 

75. P. L. Pastro, Orthogonal polynomials and some q-beta integrals of Ramanujan, J. Math. 
Anal. Appl. 112 (1985), 517-540. 

76. F. Peherstorfer, On Stieltjes polynomials and Gauss-Kronrod quadrature, Math. Comp. 55 
(1990), 649-664. 

77. , On the asymptotic behaviour of functions of the second kind and Stieltjes polyno- 
mials and on Gauss-Kronrod quadrature formulas, J. Approx. Theory (to appear). 

78. O. Perron, Die Lehre von den Kettenbruchen, Teubner, Leipzig, 1913 (band I 1954, band II 
1957); Chelsea, New York, 1950. 

79. G. Polya, Sur un theoreme de Stieltjes, C.R. Acad. Sci. Paris 155 (1912), 767-769. 

80. K. Posse, Sur les quadratures, Nouvelles Ann. Math. (2) 14 (1875), 49-62. 

81. , Sur Quelques Applications des Fractions Continues Algebriques, Hermann, Paris, 

1886. 

82. M. Prevost, Stieltjes type and Geronimus type polynomials, J. Comput. Appl. Math. 21 



IMPACT OF STIELTJES' WORK 



29 



83. A. Pringsheim, Uber die Konvergenz unendlicher Kettenbruche, S.B. Bayer. Akad. Wiss. 
Math.-Nat. Kl. 28 (1898), 295-324. 

84. , Uber einige Konvergenzkriterien fur Kettenbruche mit komplexen Gliedern, S.B. 

Bayer. Akad. Wiss. Math.-Nat. Kl. 35 (1905), 359-380. 

85. M. Riesz, Sur le probleme des moments, Ark. Mat. Ast. Fys. 16 (1921), no. 12 (23 pp.); 16 
(1922), no. 19 (21 pp.); 17 (1923), no. 16 (52 pp.). 

86. , Sur le probleme des moments et le theoreme de Parseval correspondant, Acta Litt. 

Acad. Sci. (Szeged) 1 (1922-1923), 209-225. 

87. L. J. Rogers, Second memoir on the expansion of certain infinite products, Proc. London 
Math. Soc. 25 (1894), 318-343. 

88. , On the representation of certain asymptotic series as convergent continued fractions, 

Proc. London Math. Soc. (2) 4 (1907), 72-89. 

89. I. Schur, Uber die Verteilung der Wurzeln bei gewissen algebraischen Gleichungen mit ganz- 
zahligen Koefficienten, Math. Z. 1 (1918), 377-402. 

90. A. Selberg, Bemerkinger om et multipelt integral, Norsk Mat. Tidsskr. 26 (1944), 71-78. 

91. J. A. Shohat, E. Hille and J. L. Walsh, A Bibliography on Orthogonal Polynomials, Bulletin of 
the National Research Council, vol. 103, National Research Council of the National Academy 
of Sciences, Washington, D.C., 1940. 

92. J. A. Shohat and J. D. Tamarkin, The Problem of Moments, Amer. Math. Soc, Providence, 
RI, 1963. 

93. C. L. Siegel, The trace of totally positive and real algebraic integers, Ann. of Math. 46 (1945), 
302-312. 

94. H. Stahl and V. Totik, General Orthogonal Polynomials, Encyclopedia of Mathematics and 
its Applications, Cambridge University Press, Cambridge, 1992. 

95. T. J. Stieltjes, Quelques recherches sur la theorie des quadratures dites mecaniques, Ann. 
Sci. Ecole Norm. Sup. (3) 1 (1884), 409-426. 



96. , Note a V occasion de la reclamation de M. Markoff, Ann. Sci. Ecole Norm. Sup. (3) 

2 (1885), 183-184. 

97. , Sur quelques theoremes d'algebre, OR. Acad. Sci. Paris 100 (1885), 439-440. 

98. , Sur les polynomes de Jacobi, OR. Acad. Sci. Paris 100 (1885), 620-622. 

99. , Sur certaines polynomes qui verifient une equation differ entielle lineair du second 

ordre et sur la theorie des fonctions de Lame, Acta Math. 6 (1885), 321-326. 

100. , Sur les racines de I'equation X n = 0, Acta Math. 9 (1886), 385-400. 

101. , Sur la reduction en fraction continue d'une serie procedant suivant les puissances 

descendantes d'une variable, Ann. Fac. Sci. Toulouse 3 (1889), Hl-17. 

102. , Sur la valeur asymptotique des polynomes de Legendre, OR. Acad. Sci. Paris 110 

(1890), 1026-1027. 

103. , Sur les polynomes de Legendre, Ann. Fac. Sci. Toulouse 5 (1890), Gl-17. 

104. , Sur les racines de la fonction spherique de seconde espece, Ann. Fac. Sci. Toulouse 

4 (1890), Jl-10. 

105. , Recherches sur les fractions continues, Ann. Fac. Sci. Toulouse 8 (1894), Jl-122; 9 

(1895), Al-47. 

106. , Sur certaines inegalites dues a M. P. Tchebychef (Euvres Completes de Thomas 



Jan Stieltjes, vol. II [2] , pp. 586-593. 

107. M. H. Stone, Linear Transformations in Hilbert Space and their Applications to Analysis, 
Amer. Math. Soc. Colloq. Publ., vol. 15, Amer. Math. Soc, New York, 1932. 

108. G. Szego, Bemerkungen zu einer Arbeit von Herrn M. Fekete: "Uber die Verteilung der 
Wurzeln bei gewissen algebraischen Gleichungen mit ganzzahligen Koefficienten" , Math. Z. 
21 (1924), 203-208. 

109. , Uber gewisse orthogonale Polynome, die zu einer oszillierenden Belegungsfunktion 

gehoren, Math. Ann. 110 (1934), 501-513. 

110. , Orthogonal Polynomials, Amer. Math. Soc Colloq. Publ., vol. 23, Amer. Math. 

Soc, Providence, RI, 1975 (4th edition). 

111. M. Tsuji, Potential Theory in Modern Function Theory, Chelsea, New York, 1959. 

112. J. L. Ullman, On the regular behaviour of orthogonal polynomials, Proc. London Math. Soc. 
(3) 24 (1972), 119-148. 

113. W. Van Assche, Orthogonal polynomials, associated polynomials and functions of the second 



30 



WALTER VAN ASSCHE 



114. E. B. Van Vleck, On the polynomials of Stieltjes, Bull. Amer. Math. Soc. (2) 4 (1898), 
426-438. 

115. , On the convergence of continued fractions with complex coefficients, Trans. Amer. 

Math. Soc. 2 (1901), 215-233. 

116. , On an extension of the 1894 memoir of Stieltjes, Trans. Amer. Math. Soc. 4 (1903), 

297-332. 

117. N. Ya. Vilenkin, Special Functions and the Theory of Group Representations, Translations 
of mathematical monographs, vol. 22, Amer. Math. Soc, Providence, RJ, 1968. 

118. G. Vitali, Sopra le serie di funzioni analitiche, Rend. R. 1st. Lombardo 36 (1903), 771-774; 
Ann. Mat. Pura Appl. 10 (1904), 73. 

119. H. S. Wall, Analytic Theory of Continued Fractions, Chelsea, Bronx, NY, 1973. 

120. H. Widom, Polynomials associated with measures in the complex plane, J. Math. Mech. 16 
(1967), 997-1013. 

121. , Extremal polynomials associated with a system of curves in the complex plane, Adv. 

in Math. 3 (1969), 127-232. 

122. S. Wigert, Sur les polynomes orthogonaux et I 'approximation des fonctions continues, Ark. 
Mat. Astron. Fysik 17 (1923), no 18 (15 pp.). 

123. H. Zaheer, Stieltjes and Van Vleck polynomials, Proc. Amer. Math. Soc. 60 (1976), 169-174. 

Department of Mathematics, Katholieke Universiteit Leuven, Celestijnenlaan 
200 B, B-3001 Leuven (Heverlee), BELGIUM 
E-mail address: fgaee03@ccl.kuleuven.ac.be 



